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MULTI-RESOLUTION ANALYSIS FOR ENO SCHEMES 


Ami Harten 1 

School of Mathematical Sciences, Tel-Aviv University 

and 

Department of Mathematics, UCLA 

ABSTRACT 

Given a function u(x) which is represented by its cell-averages in cells which are formed 
by some unstructured grid, we show how to decompose the function into various scales of 
variation. This is done by considering a set of nested grids in which the given grid is the 
finest, and identifying in each locality the coarsest grid in the set from which u(x) can be 
recovered to a prescribed accuracy. 

We apply this multi-resolution analysis to ENO schemes in order to reduce the number of 
numerical flux computations which is needed in order to advance the solution by one time- 
step. This is accomplished by decomposing the numerical solution at the beginning of each 
time-step into levels of resolution, and performing the computation in each locality at the 
appropriate coarser grid. We present an efficient algorithm for implementing this program 
in the one-dimensional case; this algorithm can be extended to the multi-dimensional case 

with cartesian grids. 
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1. Introduction 

In this paper we consider the Initial-Boundary Value Problem (IBVP) for hyperbolic 
systems of conservation laws in s-space dimensions. 

+ div/(u) = 0, xeVcM 3 , t>0 ( Lla ) 

u(x,0) = u 0 (x), xeV (^Ab) 

with given boundary condition on dV, the boundary of V. We assume that the problem is 
well-posed and denote its evolution operator by E(t)\ note that it includes the influence of 
the boundary conditions. 

The computational domain T) is divided into cells { Cj } 

V = U jCj, Cj r\Ck = <f> for j ^ k ( 1,2a ) 


and we assume that there is a refinement parameter h such that the largest sphere contained 
in each of the cells is of radius 0(h) , and that the ratio between the largest cell to the 
smallest one in the computational domain remains bounded under refinement. 


Let \Cj\ denote 

M-Jo* 

and let Uj denote the cell-average of u(x) over Cj 


(1.26) 


u 


j 



A(Cj) ■ u(s); 


(1.3) 


here A(Cj) denotes the cell- aver aging operator. 

Given cell-averages u = {*,) of u(x) in ©, we denote by an r-th order piecewise- 

polynomial reconstruction of u from u, i.e., 

R(x]u) = Ri(x-,u) for xtCi (1.4a) 

where Ri(x)u) is a polynomial of degree r - 1. Expressing Ri(x]u) as a finite Taylor series 
around the centroid c,- = A(Ci) ■ x 

jj j ( I1 «) = ©o + Ep£(*-VA 

fc=i K • K|=fc 


where 
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'‘ |,|£>, = '‘ W |? (c ' ) + 0(( ‘ r) ’ (1.4c) 

(accuracy) 

D 0 = u { - J2 Tf [ A (Ci) • (ar - c,-)*] D t (1.4c?) 

k=i K - |/|=* 3 


(conservation). 

Note that (1.4d) implies that A(Ci) ■ Ri(x;u) = In (1.4) we have used a multi-index 
notation 

*=(*!, •”,*.)» W=/l +•••"+/, (4>0) 
with the standard convention 


/ = (yiY 1 ~ - 


aw 


dx* dxi • • • dx *• 

We consider the numerical solution of (1.1) by the class of schemes 


Vj +1 = A(C 3 )E{t)R ( •; v n ), v] = A(Cj)u o(x) 


(1.5a) 


where v 1 - is an r-th order approximation to the cell-average of the solution u at time t n 


v j ~ A {Cj) -u(x,t n ). (1.56) 

Due to the divergence-free form of the PDE (1.1), the scheme (1.5a) takes the conservation 
form 

” r ' mVl -AJoL, mt) - R( --' v ^- Nisit (IJfe) 

where dCj is the boundary of the cell Cj and N is its outward normal. We refer the reader 
to [3] for details. 

The purpose of this paper is to present some preliminary results regarding the application 
of multi-resolution analysis to the numerical solution of hyperbolic systems of conservation 
laws. Typically these solutions contain discontinuities (shocks, sliplines, material boundaries, 
combustion fronts) which may move around and also some localized high-frequency smooth 
behavior which is associated with shedding of vortices. In such situations we have to use a 
fine grid in order to resolve the details of discontinuities and vortices, while in other parts of 
the solution the use of a coarser grid is adequate. 

The traditional solution to this computational problem is to use a nonuniform adaptive 
grid where at each time-step the discretization points are redistributed with the goal of min- 
imizing the truncation error of the scheme. In the context of the schemes (1.5) this means 
that the cells Cj = Cj 1 are redefined at each time and that their size is highly nonuniform: 
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they are very small around discontinuities and rapid smooth variation of the solution, and a 
lot larger elsewhere. This is certainly a reasonable solution to this computational problem 
and there are several such computer codes that accomplish this goal successfully. How- 
ever, the adaptive grid approach is complicated and the programming effort reqmred for it. 

implementation is formidable. . . , 

In the following we describe our concept of multi-resolution analysis which is borrowe 

from the realm of image-compression techniques. There the functions are defined on a 
uniform fine grid and are assumed to be over-resolved in some parts of it. The purpose 
of the multi-resolution analysis in image compression is to determine appropriate levels of 
resolution for the various parts of the image in order to eliminate superfulous information. 
In the context of the numerical solution (1.5) this means that our computational grid (1.2) is 
rather uniform with respect to the parameter h and is assumed to be fine enough to capture 
all the details that we are interested in, and that the solution is in fact over reso ve on 
this fine grid in large parts of the computational domain. In order to apply multi-resolution 

analysis we construct an hierarchy of nested grids as follows: 

Given cells [C]} of size h k and cell-averages (1.3) of u(x) on them, we define a 

coarser grid {C k+l } of size h k+1 > h k by joining some of the cells in the k - th grid into a 

single larger cell . 

c fc+i =u . C fc (l-6a) 


and define 




( 1 . 66 ) 


The original grid is thus the finest in the hierarchy, and we denote its cells by {C°} and their 
size by h 0 = h. To each of these grids we can apply the reconstruction (1.4) which we denote 
by R k {x;u) and also apply the numerical scheme (1.5) in order to advance its solution in 

tim What we mean by multi-resolution analysis is the assignment k(i) such that 

\R k{i \x\u) - R°(x-,u)\ < £ for xeCf ( L7 ) 

where fc(i) is the largest integer for which (1.7) holds. In Section 2 we shall describe a method 

to determine k(i). , 

The general idea is to use a multi-resolution reconstruction R(x,u ) w ic is e ne on 

the finest grid by 

R(x-,u) = for xeC? V**) 

and to consider the numerical scheme 


u? +1 = A{C°)E(T)R(-,v n ). 


( 1 . 86 ) 
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We observe that the scheme (1.8) can be thought of as an adaptive-grid calculation in which 
the cells are a combination of cells of various scales. We would like to avoid this interpretation 
because of its inherent complexity, and we look for a simple implementation. 

In Section 3 we develop such an implementation by studying a related problem and in 
Section 4 we describe an efficient multi-resolution ENO scheme for the one-dimensional case. 

In Section 5 we summarize the present state of affairs and point out directions for future 
development. 

2. Determination of Resolution Levels 

In this section we describe how to decompose a given set of cell- averages {iZ,} of a function 
u(x) over cells { Cj } of size h into a set of scales 

h = h 0 < hi < • • • < fix (2.1a) 

given by the hierarchy ( 1 . 6 ). This is done by determining the largest k = k(i ) so that 

\R k {x',u) — i?,(a:;u)f < e for xeCi (2.16) 

for a prescribed level of accuracy e. Here R,(x; u) is the polynomial of degree r - 1 in (1.4) 
and R k (x-,u) is the reconstruction which is associated with the scale h k . 

First let us describe the way in which the reconstruction R(x;u) is defined for the given 
grid. Expressing Ri(x\u ) as a finite Taylor series around the centroid c, (1.4b) we have to 
find coefficients {Tb},0 < \£\ < r - 1 which satisfy the requirements (1.4c) - (1.4d). Let us 
denote by d the vector of unknowns 

(2.2) 

which are ordered in groups of equal \£\ with \£\ = 0, 1, • • • , r - 1, and denote by « = «(r) the 

number of unknowns in (2.2). In [3] we show how to select a stencil of k cells with indices 
J(i) (including i ) so that the relations 

^(Cy) • Ri(x] u ) = uj, j e J(i) (2.3 a) 

result in an invertible system of linear equations 


Here Q is the matrix with entries 


Qd = u. 


(2.36) 


A(Cj) • (x - cif/h* 


(2.4a) 
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and u is an appropriate ordering of the RHS of (2.3a). Note that due to the scaling in (2.2), 
the entries of the matrix Q (2.4a) are 0(1) under refinement and consequently 

||Q _1 II — cons ^ as h 0. (2.46) 


In [3] we consider two cases: (1) a fixed choice of a centered stencil J(t) winch results in 
an upwind biased “linear” scheme (1.5); (2) an adaptive choice of stencil J(.) - 'W® ) 
where the stencil selected is the one in which u(x) is the smoothest among several candidate 
stencils. This adaptive selection results in an ENO scheme (1.5). 

We turn now to describe the decomposition into scales (2.1). For k = K, K - 1, 

check whether 

max | A[Cj) • R k {x’, u) — Uj\ < 6(e); 


we 


(2.5) 


Mi) is the largest k for which (2.5) holds. Next we show that (2.1b) holds for an appropriate 

choice of 6(e). Let us rewrite the polynomial of degree (r - 1) R l (x l«) «s a f ' nl,e Tay or 
series around the centroid « as in (1.4b) and denote the vector ordering of its coefficients 

(2.2) by d. We denote by u k the vector ordering of 



u k = A{Cj) • R k (x-,u), jtJ{i ) 

(2.6a) 

as in (2.3a). Clearly 

Qd = u k 

(2.66) 

where Q is the matrix in (2.3b). It follows therefore that 



Q(d -d) = u-u k , 

(2.6a) 


|| d - d\\ < HQ' 1 II ||« - u k \\ < ||Q” 1 ||6(e) 

(2.7a) 

and consequently _ 


max \R k (x\u) - Ri(x]u ) | < C'||Q- 1 ||6(e) 

xtC{ 

v z - ' 

where 


(2.7c) 

\x c t \ 

C = max max , ... 

i<KI< r_1 xtCi 6 1 * 1 

and || |l is the £ r norm. 

We note that the terms C and || Q"MI ( 2 - 7 ) de P end on 

the geometry 


of the cell and of the stencil J (i), respectively. 

Up to this point we have dealt with the concept of multi-resolution analysis in the context 
of approximation of functions. We turn now to consider its application to the numerical 
solution of the IBVP (1.1) by the scheme (1.5). Given {«*}, the cell-averages of the numerical 
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solution at time t n over the cells {Cj}, we apply the multi-resolution analysis (2.5) to find 
k(i) and define the multi-resolution reconstruction R(x;v n ) (1.8a) by 

R(x;v n ) = R k ^(x-,v n ) for xe (7, ( 2 .8 a ) 

and consider the numerical scheme 


vf' = A{Ci)E{T)R{x-,v”). 

We observe that due to the well-posedness of the problem 


( 2 . 86 ) 


||£(r)f?(x; v n ) - E(t)R{x ; u n )|| < const. ||R(x; v n ) - R ( x ; u n )|[ (2.8c) 

and therefore the deviation of the values computed by the multi-resolution scheme (2.8) 
from those of the original scheme (1.5) is of order e. Roughly speaking, all the values vj» 
m (2.8b) for which k(j) = k 0 can be obtained from the calculation on the coarse grid k 0 by 

V F' = ^(C,)£(r)R*°(x;u n ), k(j ) = k 0 . ( 2 .9) 

What we need now is an efficient algorithm that will enable us to perform the computation 
in (2.8) using the appropriate coarse grid calculations (2.9) at the computational cost of a 
corresponding adaptive grid implementation without its inherent complexity. To get ideas 
for such an algorithm we consider a related computational problem in the next section. 


3, A Related Problem 


In this section we consider the problem of computing discrete values of a composite 
function (f>(v(x)), where v(x) is a function of the type that we encounter in the solution of 
hyperbolic conservation laws (1.1), he., it is piecewise smooth with different scales of smooth 
variation. The function <j>{v) is a model for the numerical flux in (1.5c). We assume that 
is a smooth function of v which is expensive to compute, e.g., 

/ OO 

g(v,x)dx ( 3 .i) 

with g(v,x) that depends smoothly on v. 

In this model problem we assume that u(x) is defined in [0,1] and its discrete values are 
given on the uniform grid 

G° = { x i))=Q, Xj=2~ J -j, Vj = v(xj). (3. 2a ) 

We assume that the grid G° is fine enough to resolve v(x) to our satisfaction, and that in 
fact u(x) is over-resolved in some parts of the grid. Our task is to find an efficient algorithm 
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2 J which approximate ^(v(*j)) within a specified tolerance 


(3.26) 


to obtain values 0 < j < 

of error e, i.e., 

| $ - <i>(vj)\ < e- 

We consider the nested set of grids {0*}, 0<k< AT, which is defined by 

g‘ = {* j -w, J k = Vi)£o (3 - 3a) 

in which G° is the finest (given) grid and G K is the coarsest; let h k denote the spacing of 
points in the grid G k 


h k = 2 


,k-J 


(3.36) 


Clearly in 

G k n G fc_1 = G k = {x 2 /}2 itj*-' > * ' ’ 

G k ~ x -G k - {x 2 <-i}(2<-i)eJ*- 1 - (3.3d) 

The restriction v t = of the values {v,} to the grid & constitutes a level of 

resolution of v(x) which corresponds to the scale h k . 

We begin the computation on the coarsest grid G K by calculating the values } jiJ K 

from the given “expensive” expression, say (3.1) 

df = «»,•). j<J K - (3A) 

Then for k = K 1 we use the already computed values on the grid G* to obtain 

acceptable approximate values **-> of *(»(*)) »■> «•<* * rid G ‘~' ' " f f° WS; 

At all points (3.3c) which are common to both grids, we retain the already compute 

^ #*-*. 2W~. < 35 °> 

The points (3.3d) are the ones added to G‘ by splitting each of its intervals into two. At 
each such middle point we make a decision whether to compute 4> directly, i.e., 

= «»*-.), (2f - (3 56) 

or to approximate (f> there by interpolation from the grid G l , i.e., 

fah = /(*«-« /), (2£ — l)ej fc_1 . 

We elect to interpolate (3.5c) wherever it can be ensured that the interpolated value of 4> is 
accurate to a prescribed tolerance e k -i 

(3.5 d) 


(3.5c) 


— <K u 2*-l)l < £k-U fit 


fc-1 
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The basic idea behind this algorithm is that is to be computed directly until the 
prescribed level of resolution implied by (3.2b) is obtained. Once the required level of 
resolution is achieved in a certain locality we cease to use the expensive direct computation 
and the values of <S are transferred to finer grids by interpolation (which is assumed to be 
consi erably less expensive). The basic assumption is that 4>(v) is smooth and does not 
vary muc in [mmv(x),maxv(x)]. Thus the quality of approximating Mvlxl) on G*“ 1 by 
interpolation from the coarser G* is determined up to a scaling factor by the error in doing 
so or v(x). Since v(x) is given on all grids, our success in interpolating «*-* from u* can be 
easily measured. The criterion for acceptable interpolation is formulated by 

max|/(^ + «ti‘)-„ i+ ,| \I(x,\<t>(v k )) - <f,{v,)\ <e k (3.6a) 

Here the “max" denotes checking the quality of the interpolation at the point jej k -' (t = 0) 
and possibly at neighboring points of the finest grid (/ = ±1, . . .). St is the corresponding 
tolerance for the u-interpolation, and it is to be determined by analytical methods. 

We recall that the values of ** that are actually used in the interpolation (3.5c) are an 
£k approximation to the exact values <f>(v k ), thus 


\I(x-,4> k ) I(x] <j>(v k ))\ < Ck • £k, ( 3 . 66 ) 

and therefore in (3.5d) 

|/(x j; **) - < |/(*,; 4> k ) - J( Xj] *(»*))| + I /(*,.; ^(«*)) - # VJ )I < (1 + c„) ■ e h . (3.6c) 

Consequently the tolerance levels {e*} satisfy 


£o - e > ei > • • • > e K = 0, (3.7a) 

(1 + Ck) ■ £k < £k-i- ( 3 . 76 ) 

We turn now to consider the practical implementation of the multi-resolution algorithm 
(3.4) - (3.6). The computational effort of executing it is composed of fV(l - 2~ K ~ 1 ) checks, 
M < N direct calculations of <f> and (N - M) interpolations of <j>; N = 2 J . The number 
of checks is fixed and is almost N, so we want to make the test in (3.6a) as simple as 
possible. On the other hand we wish to reduce M, the number of direct evaluations <t>(v)- 
this calls for an elaborate test so as to not miss points in which interpolation is acceptable.’ 
hese considerations lead to a compromise in the level of sophistication to be used in the 
interpolation I(x,v ) and J(x; <f> k ). The simplest and least expensive interpolation is based 
on the use of a fixed central stencil. However, v(x) is a discontinuous function and using 
a fixed central stencil means that we shall miss all the points near the discontinuity which 


otherwise could be well-approximated by the more expensive ENO interpolation that uses 
an adoptive stencil [1], [2]. Investing even more in an ENO interpolation by using subcell 
resolution [4] will enable us to get a good approximation even in a cell which contains a 
discontinuity. 

A reasonable strategy is to use a hierarchy of checks: first to try the simplest interpolation 
with fixed central stencil which is a good enough test for most of the domain. If the test in 
(3.6a) fails for the simple interpolation we may try to get an acceptable approximation by an 
ENO interpolation. In cases where the computation of <j>(v) is so expensive that it justifies 
an additional investment in checks, we may try to use also a subcell resolution approach. 

4. Multi-Resolution for One-Dimensional Conservation Laws 

In this section we present a multi-resolution algorithm for the solution of one-dimensional 
hyperbolic conservation laws which is motivated by the algorithm (3.4) - (3.6) for the model 

problem. 

We consider the IB VP for the one-dimensional conservation law 

u t + f(u)* = 0, 0 < x < 1, t> 0 (4-la) 

u(x,0) = ti 0 (x), 0<x<l (4-1&) 

with appropriate boundary conditions at x = 0,x = 1. 

We discretize functions in [0,1] by taking their cell-averages on the grid G° (3.2a) 

Tij = A([xj_i,Xj]) • u(t), 1 < 3 < 2 J - 

Given cell-averages H = {uj}j c j° we consider the reconstruction R(x] u) which is defined 


in a piecewise manner in [0,1] and satisfies 

R(x-, u) = u(x) + 0 (h r ), wherever u(x) is smooth, (4.2a) 

^([xj-i, Xj]) • R(x;u) = Uj, 1 <j< 2 J - ( 4 - 26 ) 

We use the numerical scheme (1.5) which can be written in this case as 
vf 1 = Adxj^xj}) • E{t) • ii(x; v n ), v'j = ^([x^.x,-]) • «o(*), 1<J< 2 J . ( 4 -3a) 

and has the conservation form 

”J +1 = - ^o(7j - 7j-i)> >‘o = r/h 0 (4.36) 

with the numerical flux 

7, = - [' /(«(<) ■ (4 3c) 

1 T JO 


9 


Note that the computation in (4.3) is to be performed with A 0 which is restricted by the 

CFL condition, and that f 0 and f 2 j include the influence of the corresponding boundary 
conditions at x = 0 and x = 1. 

In solving one-dimensional problems, as well as in multi-dimensional problems on carte- 
sian grids, it is convenient to work with the primitive function 


U{x,t) = j u(£,t)dt. ( 4>4a ) | 

Since 1 

-^[U(xj,t n ) — U(xj-i,t n )] = vj , (4.46) Z 

~[U{x 3 , t n + T )~ U( Xj , t n )] = (4.4c) I 

it is convenient to define the reconstruction R(x } v n ) by i 

m 

R(X ’ = li 1 ^' U ") ( 4 -5) I 


where /(*; £/") is a piecewise interpolation of {{/(*„<„)} (see [2]) and to compute the nu- 
merical fluxes from (4.4c). 

We turn now to describe a multi-resolution algorithm for carrying out the computation 
of the numerical solution (4.3). We are given v n = {v?} jtJ0 on the grid G° and our task 
is to compute all the numerical fluxes {/_,} on this grid within a prescribed tolerance for 
error. Our basic assumption is that the solution is well-resolved on the grid G°, and in fact it 
may be over- resolved in some parts of it. The computational task is thus similar to the one 
entailed in the model problem, except that here J j is really a functional of u(x, t n ) and not a 
function like (f>{v(x )) in Section 3. We use the same set of nested grids (3.3) G k , 0 < k < K, 

and initialize the algorithm by evaluating ~J ■ in (4.3c) on the coarsest grid G K . These fluxes, 
which we denote by 

fj = 7K)j, jeJ K , (4.6) 

are actually an r-th order approximation to (4.3c) which uses the reconstruction f?(x;u n ) of 
the finest grid G°. Thus we use the same numerical fluxes that we would use in a straight 
forward fine grid calculation, except that these fluxes are computed only in the few locations 
which correspond to the coarsest grid G K . 

The algorithm for the numerical flux computation proceeds in analogous way to (3.5). 
For k = K, . . . , 1 we use the already computed values of the numerical flux J k on the grid 
G to obtain acceptable values f k 1 on the finer grid G k ~ l as follows: 

At all points (3.3c) which are common to both grids we retain the already computed 
values 

7‘;‘ = 7*<> 2&J*-*. (4.7a) 
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At the new middle points (3.3d) we make a decision whether to compute the numerical 
flux directly . 

f-\ = /(»*)»_„ (M-1W*- 1 ( 4 - 76 ) 

or to approximate it by interpolation of f, the already computed values of the numerical 
flux on the grid G k 

?£. = /(*«-. ,T). (“- (47c) 

The decision whether to interpolate or compute at a certain locality depends on whether 
we can obtain the prescribed accuracy for the numerical flux from the level of resolution 
which is offered by the grid G k . As is indicated by (2.8c), the error 

f-7, = - [’(f(E(t) • - f(E(t) ■ R\„)]dt 3 C f (4.8a) 

3 3 T JO 

which is committed by replacing the exact numerical flux (4.3c) by the numerical flux 
corresponding to the level of resolution which is available in G k , can be bounded by the 
corresponding error in reconstruction 


max | R k (x; v n ) — R{x ; u n )| — ef ■ 


(4.86) 


In the constant coefficient case f(u) = au, we get 


i £ /| < M f |i? fc ( Xj - at; v n ) - R{xj - at ; v n )\dt < \a\ef, (4.9) 

and similar estimates can be obtained in the general nonlinear case. As is indicated by (2.5) 
and (2.7), ef (4.8b) can be estimated from e" 

e u = max |A([*i_i,*i]) ■ H*(*;w n ) - w"l* (4.10a) 

3 . £ J(i)uJ(i-i) 

This check is foolproof but expensive. In our numerical computations we have experimented 
with a simplified version of (4.10a) 

= max{|A([z j -„*i]) • -»?|, |X([*i,*i+.]) ■ $(* ;»’) “ »’«l) t 4 ' 106 ) 

which takes into account only the two adjacent fine grid intervals that determine the numer- 
ical flux at Xj. These numerical experiments indicate that the test (4.10b) may be adequate 

enough. 

From our numerical experiments we have also learned that it is better not to interpolate 
the numerical fluxes directly (4.7c), but rather use the following indirect way: Given / we 
can compute v n+1 ' k , the cell-averages at time f„+i on G k by 

v] +hk = vf - ^(7- - 75-2*). ( 4 - lla ) 
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From v n+1 ’* we get the values of the primitive function U n+1,k by using the relation (4.4b). 
Using interpolation 7(x; U n+1 ’ k ) for the primitive function at time f n+1 and possibly a dif- 
ferent interpolation 7(x; U n,k ) for the primitive function at time t n we can now obtain the 
value of the numerical flux at the middle points (3.3d) from the relation (4.4c). Thus we 
replace (4.7c) by 

TjT-1 = “[/(**- 15 U n+1 ’ k ) - U n ’% (2£ - 1 )cJ k -\ (4.116) 

We remark that unlike the model problem of Section 3, the computation of the numerical 
flux involves propagation of discontinuities. Consequently a stencil of cells which did not 
contain discontinuities of u(x, t n ) may contain a discontinuity of u(x, t n+1 ) at its two extreme 
cells. To overcome this problem we use a central stencil for R k in the test (4.10b) and for the 
interpolation 7 for the primitive function at time < n ; in order to account for the movement 

of discontinuities we take 7 in (4.11b) to be an ENO interpolation for the primitive function 
at time 2 n+1 . 

5. Summary 

This work on multi-resolution analysis has been motivated by the current interest in 
wavelets. The concept of wavelets has developed into a beautiful theory which is rich in 
structure. The endowment of wavelets with many properties is also their main handicap when 
it comes to practical computations: Indeed wavelets are more local than Fourier analysis, 
but y are net really local. The desire to have an orthonormal basis with wavelet functions 
of compact support results in wavelet functions which are of fractal nature. The attempt to 
make these wavelet functions smoother results in considerably enlarging their support. 

In the present work we have attempted to take from wavelet only its most important 
idea - the multi-resolution analysis - and to strip it of any other structure. In doing so we 
have developed a multi-resolution analysis which is devoid of beautiful theory, but we have 
gained flexibility and the ability to adapt locally. Given a function u(x) which is represented 
by cell-averages on an arbitrary partition into cells, we have shown how to decompose it into 
various scales. This general method provides a very useful tool for removing superfulous 
information, especially when the data is only piecewise smooth, e.g., application of this 
methodology to image-compression has proven to be very successful. 

We remark that this multi-resolution analysis can be easily extended to functions u(x) 
which are represented by point-values on an unstructured grid. The reason for presenting 
only the cell-average version is that we wish to combine multi-resolution analysis with high 
order of accuracy, and this forces us to choose the scheme (1.5) for accurate cell-averages. 
We cannot use point-value schemes which are both high-order accurate and in conservation 
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form. The need to fictitiously write f x in (4.1a) as a difference of a numerical flux results in 
a strong dependence of the numerical flux of the point-value scheme on the spacing of the 
grid. Consequently we cannot retain the already computed numerical fluxes on the grid G 
when refining to G*' 1 (4.7a) for point-value schemes with order of accuracy larger than 2. 

The multi-resolution analysis presented in this paper has points in common with both 
multi-grid and adaptive-grid methods. It resembles multi-grid in its use of nested grids. 
However multi-grid is primarily an iteration scheme which is designed to accelerate con- 
vergence, while multi-resolution is primarily a data-compression scheme which is designed 
to avoid superfulous computations by discarding insignificant information. Its goal is the 
same as that of adaptive-grid methods but it differs in its concept: Adaptive-grid methods 
tell you where to refine the mesh, while multi-resolution analysis tells you where not. The 
difference in the practical implementation of the two concepts reminds one of the difference 
between shock-fitting and shock-capturing. In adaptive-grid methods you have to worry 
about missing spontaneously generated features. On the other hand you can do an optimal 
job in resolving features that you know of, although at the cost of complex programming. 
In multi-resolution methods you cannot miss new features because conceptually we always 
compute on a fine enough grid. The main issue is efficiency: You want to make sure that 
you do not compute on too fine a grid, when you do not have to. An additional advantage is 
the automatic way in which the data is handled in multi-resolution analysis which translates 
into simplicity of programming. 

The one-dimensional algorithm of Section 4 can be extended in a straightforward way to 
the multi-dimensional case with cartesian grids. Based on our one-dimensional results we 
are led to believe that realistic savings by factors between 10 to 100 are quite possible in the 
three-dimensional case. 
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